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By numerically integrating the three-dimensional Maxwell equations in the time domain with reference to a 
dispersive quadratically nonlinear material, we study second harmonic generation in planar photonic crystal 
microrcsonators. The proposed scheme allows efficient coupling of the pump radiation to the defect resonant 
mode. The out-coupled generated second harmonic is maximized by impedance matching the photonic crystal 
cavity to the output waveguide. 

In the early sixties the very first experiments in nonlin- 
ear optics were Frankens investigations of traveling- wave 
second harmonic generation. 1 It took a number of years 
to realize that guided-wave optics and resonators could 
dramatically enhance phenomena requiring high inten- 
sities and tuning of critical parameters. 2-4 Since then, 
both the advent of nano-optics and the technological ad- 
vances have made new solutions available. Among them, 
microresonators can certainly be considered among the 
best candidates for nonlinear optics and frequency gen- 
erations. 5-8 In the past few years, photonic crystal (PC) 
microcavities, i. e. periodic (bandgap) structures host- 
ing a resonant defect, have attracted attention for some 
of their unique properties. 9, 10 In particular, in PC mi- 
croresonators it is possible to obtain extremely high qual- 
ity factors (Q) in reduced volumes while tailoring their 
dispersive features, 11, 12 characteristics which can be ex- 
ploited to achieve efficient frequency generation with var- 
ious schemes. 13-16 The advantages inherent to a high 
confinement, however, are partially counterbalanced by 
a difficult energy coupling into the resonators: well iso- 
lated resonant states, in fact, correspond to large (exter- 
nal) quality factors of the cavity. 17 

To this extent, impedance matching has been proposed 
based on a properly designed coupling to/from the mi- 
crocavities. 18, 19 In this Letter we propose and investigate 
an efficient outsourcing scheme to maximize frequency 
doubling from a PC defect. Resorting to a second-order 
nonlincarity in a large Q photonic crystal microcavity, 
an optical pump is up-converted to a resonant cross- 
polarized harmonic signal. Since the device is designed 
to be nearly transparent to the pump wavelength, input 




Fig. 1. Artist's sketch of the structure. 

coupling losses are minimized while the out-coupled sec- 
ond harmonic can be maximized by impedance match- 
ing. 

The temporal evolution of the resonant mode at the 
second-harmonic (SH) can be described and related to 
the cavity parameters by coupled mode theory in the 
time domain (CMT-TD): 17 



da Q 



dt 



(1) 



where a SH is the mode amplitude at the SH frequency 
uj sh , t sh takes into account both internal (1/t ) and 
external (l/r e ) losses (1/t sh = l/r + l/r e ), k is the 
(three-dimensional) overlap integral between the sus- 
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Fig. 2. Quality factors Q e (diamonds) and Q a (stars) 
versus spacing d (in units of r). The inset shows the 
ratio Qe/Qo versus d/r. 
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ccptibility tensor and the electric fields at fundamental 
(FF) and SH frequencies, s FF is the amplitude of the 
FF guided mode propagating through the cavity via 
input and output waveguides. The latter amplitude 
acts as a source for SH generation in the resonator. A 
pictorial sketch of the PC microcavity with input and 
output guides is shown in Fig. 1. 

Taking s FF = £, FF exp[-iu FF t] and a SH = 
A SH exp[—iuj SH t], with uj sh = 2u> FF yields the steady 
state solution of (1) 



(2) 



The out-coupled SH radiation is linked to the external 
losses by: 



— A — 



(3) 



From the last expression, the normalized efficiency is 
then: 



2t„ 



(1 



Qo 



Qe j Qo 



(4) 

On the RHS of (4) r e and r Q have been replaced by in- 
ternal Q and external Q e quality factors of the cavity, 
respectively, being 1/Q = 2/(tu) sh ). For a constant k 
in (4), the SHG efficiency is linearly dependent on Q a 
and maximized when Q e /Qo = 1- Otherwise stated, once 
the cavity geometry is defined (i.e. given a certain Q ) 
and for small changes in the overlap integral k, the out- 
coupled SH can be optimized by engineering the factor 
Q e . In order to demonstrate the validity of the insight 
above, we need to take into account all the loss mecha- 
nisms (radiation, absorption and coupling), material and 
PC dispersion and nonlinear response. To this extent 
we employed a finite difference time domain (FDTD) 



Fig. 3. Snapshots of the electric field in the central 
portion of the PC microresonator. Top panel: FF x- 
component; bottom panel: SH y-component. 

code to study the nonlinear process without approxima- 
tions. 20,21 In the code, Maxwell equations in vacuum are 
coupled to the material polarization: 



VxE = -,i f<, VxH 
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(5) 

The last equation describes a dipole oscillation with an 
anharmonic term accounting for the material quadratic 
response. Its polarization P is due to a single-pole dis- 
persion with a Lorentzian distribution centered in cJr. 
We chose a loss coefficient 70 = 3.8 x 10 8 s _1 , ujr = 
1.1406 x 10 16 s _1 , a static permittivity e s = 11.7045 and 
uJn = 1.0995 x 10 15 s _1 to mimic the response of GaAs. 
The pertinent tensor D values were selected after a 45 
degrees rotation in the plane zy of GaAs crystallographic 
axes. This permits to generate (quasi-) TE polarized 
second harmonic when a linearly polarized (quasi-) TM 
pump-mode is launched in the structure through the in- 
put channel. We investigated a microcavity consisting of 
a PC membrane with three concentric rings of holes with 
a central defect to support a resonant state at SH. The 
radius of the PC holes is r = 90nm and the lattice con- 
stant a = 300nm. Two suspended channel waveguides 
of width w = lfim and thickness h = 500nm feed and 
out-couple the radiation to and from the microcavity, re- 
spectively The separation d can be varied in order to ad- 
just the SH impedance of the waveguide to the resonant 
mode. For the sake of simplicity, we set the separation d 
the same at both the input and the output. The structure 
is designed to be resonant at Xsh = 1.1370/xm for the 
TE polarization. The input waveguide couples a TM- 
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Fig. 4. SH output efficiency versus separation d (units of 
radius r). 

polarized mode of amplitude s FF at Xpp = 2.2740 fim 
into the resonator, where it is up-converted to the cross- 
polarized SH mode of amplitude a SH and out-coupled 
to the TE-mode s SH . As d was tuned, we verified that 
the resonant frequency did not change appreciably. With 
the FDTD code, first we linearly characterized the defect 
state. By evaluating the rate of energy decrease and the 
transmission T at CJ SH , T = (l + Q e /Qo) 2 , we calculated 
both Q and Q e for various d, according to: 



Qe = 



Qo-Q s 



h = Q 3H [i 



2(l-\/T) J 



(6) 



Fig. 2 displays the computed values: as expected, increas- 
ing the distance d between the cavity and the waveg- 
uide^) barely affects the internal (intrinsic) cavity losses, 
although an apparent increase takes place in the exter- 
nal quality factor owing to a lower coupling to and from 
the resonator. The inset of Fig. 2 shows the calculated 
ratio between Q e and Q . On the basis of Eq. Q, the 
SH efficiency is maximum when such ratio equals unity, 
i.e. for values of d/r close to one (see inset). We numeri- 
cally launched a TM-polarized cw-like excitation at ui FF 
in the feeding waveguide. The top panel of Fig. 3 is a 
snapshot of the TM component of the FF electric field, 
flowing from left to right through the PC and its defect. 
The TE-componcnt of the generated electric field at SH 
is shown in the bottom panel of Fig. 3: the generated 
mode at u) SH clearly resonates in the central portion of 
the cavity. To evaluate the SHG efficiency we computed 
the TE-polarizcd SH power in the output waveguide and 
scaled it to the square power of the FF radiation propa- 
gating through the PC. The results are displayed in Fig. 
4, which shows an apparent maximum in efficiency when 
the separation between waveguidc(s) and PC cavity re- 
alizes the impedance matching, e.g. for d/r = 1. Note- 
worthy, at variance with the simple CMT-TD, the latter 
ratio accounts for both the variations in k and the ac- 
tual FF power coupled in the structure. For this case, we 
calculated an out-coupled SH-power P SH = 0.22W when 



pumping with P FF = IKW through the input channel. 
Such value for P SH needs be linked to the small quality 
factor Q of the resonator (for computational reasons we 
analyzed a rather small PC) and to the non-optimized 
FF sourcing into the resonating mode at SH (i.e. a lim- 
ited k). Both factors could be enhanced by a proper de- 
sign. In conclusion, we have numerically investigated by 
FDTD second harmonic generation in a photonic crystal 
microcavity, maximizing the out-coupled SH power by 
a simple impedance matching based on the separation 
between the PC and the output channel. We anticipate 
that these results will promote design and realization of 
efficient integrated structures for all-optical signal han- 
dling. 
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